Enhanced spin density wave in LaOFeSb 
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We predict atomic, electronic, and magnetic structures of a hypothetical compound LaOFeSb 
by first-principles density-functional calculations. It is shown that LaOFeSb prefers a stripe-type 
antiferromagnetic phase (i.e., spin density wave (SDW) phase) to the non-magnetic (NM) phase, 
with a larger Fe spin moment and greater SDW-NM energy difference than those of LaOFeAs. The 
SDW phase is found to favor the orthorhombic structure while the tetragonal structure is more 
stable in the NM phase. In the NM-phase LaOFeSb, the electronic bandwidth near the Fermi 
energy is reduced compared with LaOFeAs, indicating smaller orbital overlap between Fe d states 
and subsequently enhanced intra-atomic exchange coupling. The calculated Fermi surface in the 
NM phase consists of three hole and two electron sheets, and shows increased nesting between two 
hole and two electron sheets compared with LaOFeAs. Monotonous changes found in our calculated 
material properties of LaOFePn (Pn=P, As, and Sb) , along with reported superconducting properties 
of doped LaOFeP and LaOFeAs, suggest that doped LaOFeSb may have a higher superconducting 
transition temperature. 

PACS numbers: 74.70.-b, 71.18.+y, 71.20.-b, 75.25,+z 



Since the discovery of a new class of superconductor 
LaOFeP with a superconducting transition temperature 
T c ~ 4 K [1] and subsequent dramatic increase to 26 K 
with LaOi-^F^FeAs compound 0], and to 43 K with 
external pressure 0], a huge amount of studies have fol- 
lowed to understand the underlying mechanism of the su- 
perconductivity in these iron-based layered oxypnictide 
compounds and to raise T c further. Often being com- 
pared with cuprate superconductors, this family of mate- 
rials is expected to give a clue to the superconductivity in 
the cuprates, with similarities such as two-dimensionality 
of structural and electronic properties, exhibition of the 
antiferromagnetic (AFM) ground state of the undoped 
parent materials, subsequent loss of the magnetism ac- 
companied by the emergence of the superconductivity 
with doping, etc. As the magnetism is believed to be 
closely related to the superconductivity of these mate- 
rials, many studies have been focused on the magnetic 
property of the materials, and a number of theoreti- 
cal studies have predicted magnetic instability in pure 
LaOFeAs @, H, @, 0, S] • It is known that there is a struc- 
tural phase transition from tetragonal to orthorhombic 
phase near 150 K 0,13, which is associated with anoma- 
lies in resistivity, dc magnetic susceptibility, and specific 
heat @,[Tl|]. At temperatures below ~ 130 K, an AFM 
order of spin density wave (SDW) type was identified 
with magnetic moments of about 0.3 [Ib @, fioj ] . The 
existence of SDW state in these materials have been con- 
firmed in many other experimental [ill . Ilil . [l3| and the- 



oretical 11, 14, 3, 16| works, and can be understood by 
a peculiar Fermi surface topology, i.e., the Fermi surface 
nesting. 

Since the conducting FeAs layer is conceived to be the 
relevant part of magnetism and superconductivity, there 
have been large efforts to tailor the structure by modi- 
fying the insulating layer, which is LaO layer in the case 



of LaOFeAs, to raise T c . The constraint considered is 
that the insulating layer should provide one electron per 
one formula unit of FeAs and substantial doping should 
be introduced by substitution or non-stoichiometry. By 
substituting La with other rare-earth (RE) elements, se- 
ries of RE(Oi_ ij ;F : j : )FeAs are reported to have T c of 55 K 
for RE=Sm 0, 41 K for RE=Ce [3, 52 K for RE=Pr 
[IH, and 50 K for RE=Nd j^. Besides the fluorine dop- 
ing, T c of 53.5 K is reported with GdOFeAs with oxygen 
deficiency [2l[. More recently, BaFe2As2, which contains 
the same FeAs layer as in the REOFeAs compounds but 
Ba layer instead of REO layer, has been found to have 
similar properties of magnetism and superconductivity, 
with the SDW transition at 140 K in undoped samples 
[2^ | . and superconducting transition at 38 K in K-doped 
ones [23|. SrFe2As2 is another similar discovery, with the 
SDW transition at 205 K [3, and superconducting 
T c of about 38 K |25J for K-doped samples. In these 
many materials, variations made on the material so far 
have been restricted to the insulating layer except for 
FeSe cases, preserving the FeAs layer itself since the first 
breakthrough of substituting P of LaOFeP with As 0. 
A straightforward idea to modify directly the FeAs layer 
is substituting As with Sb, inspired by the great success 
of As replacing P. 

In this work, we study the structural, electronic, and 
magnetic properties of a hypothetical iron-based oxypnic- 
tide compound LaOFeSb through first-principles density- 
functional pseudopotential calculations. We find that 
the SDW-type AFM spin state is the ground state for 
this compound, similar to existing other oxypnictides 
LaOFeP and LaOFeAs. The crystal structure cou- 
ples with the magnetic configuration so that the NM 
phase favors the tetragonal structure whereas the lower- 
symmetry orthorhombic structure is favored by the SDW 
state. As the distribution of Fe d bands near the Fermi 
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FIG. 1: (color online). Atomic structure of LaOFeSb. In this 
structure, LaO and FeSb layers alternate along the c direction. 
Black lines represent the tetragonal unit cell of the NM phase, 
and \f2 x v2 supercell is also indicated in gray lines. The Fe 
spin configuration in the SDW phase is shown with different 
colors of Fe atoms. 

level (Ef) is different from that of LaOFeAs, the Fermi 
surface topology is also different. The Fermi surface con- 
sists of three hole pockets around the T point and two 
electron pockets around the M points, with enhanced 
nesting between the hole and electron surfaces. We find 
that LaOFeSb in the SDW phase has a local Fe spin 
moment larger than that of LaOFeAs, making a trend 
of monotonous increase from phosphide to antimonidc, 
along with the relative stability of the magnetic ground 
state with respect to the NM solution. Our calculation 
suggests LaOFeSb, which is not yet reported to be syn- 
thesized, could be a candidate for higher T c . 

Our first-principles calculations are based on the 
density-functional theory (DFT) within the general- 
ized gradient approximation (GGA) for the exchange- 
correlation energy functional [261 ] and the ab-initio norm- 
conserving pseudopotentials as implemented in SIESTA 
code [13] ■ Semicore pseudopotentials are generated us- 
ing electronic configurations 5s 2 5p 2 5d°6s° with 5s and 
5p electrons as valence electrons, and 3s 2 3p 6 3d 6 4s° with 
3s, 3p, and 3d electrons as valence electrons, for La and 
Fe, respectively. Electronic wave functions are expanded 
with localized pseudoatomic orbitals (double zeta polar- 
ization basis set), with the cutoff energy for real space 
mesh of 400 Ry. Brillouin zone integration is performed 
by Monkhorst-Pack scheme [28[ with 12 x 12 x 6 k-point 



TABLE I: Calculated structural parameters, DOS at Ef 
(N(Ef)), and Fe magnetic moment (m) of LaOFePn, for 
Pn=P, As, and Sb. Optimized crystal structure for NM phase 
is tetragonal with 8 atoms in a unit cell and that for the SDW 
phase is orthorhombic with 16 atoms in a unit cell. Spin or- 
dering in the SDW phase doubles the unit cell size along the 
c axis with 32 atoms in a unit cell. The lattice parameters 
a and b in the orthorhombic structure are slightly greater 
than v2 times of a in the tetragonal structure. In the SDW 
phase, the Fe atoms of the same spin are located along the 
b axis. Iron (Oxygen) atoms are located at z = 0.5 (z = 0) 
along the c-axis. 







NM (Tetragonal, P4/nmm 






Pn 


a (A) 


6(A) 


c(A) 


2 (La) 


z(Pn) 


N(E f ) 
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3.972 


3.972 


8.590 


0.153 


0.624 


4.5 


As 


3.999 


3.999 


8.706 


0.145 


0.640 


1.7 


Sb 


4.106 


4.106 


9.311 


0.130 


0.659 


2.9 






SDW (Orthorhombic, Cmma) 




Pn 


a D (A) 


b (A) 


c(A) 


z (La) 


z(Pn) 


m(p B ) 


P 


5.741 


5.672 


8.641 


0.150 


0.634 


2.42 


As 


5.780 


5.693 


8.875 


0.139 


0.654 


2.83 


Sb 


5.955 


5.844 


9.542 


0.124 


0.673 


3.13 



grid. The unit-cell parameters are optimized by the total 
energy minimization and used for the electronic structure 
calculations. Atomic forces are relaxed below 0.04 eV/A. 
When our calculational method is applied to LaOZnPn 
(Pn=P, As, and Sb), as a test, the obtained structural 
parameters are in good agreement with experimentally 
known values. 

Oxypnictide compounds have layered structure with 
alternating LaO and FePn layers along the c-axis, as 
shown in Fig. 1. In the FePn layer Fe atoms form a 
two-dimensional square lattice and are four-fold coordi- 
nated with Pn atoms to form a network of tetrahedrons, 
and the LaO layer is in structural analogy with FePn 
layer with Fe and Pn atoms replaced by O and La atoms, 
respectively. LaOFeAs is known to have two structural 
phases of different s ym metries with a transition temper- 
ature of ~ 150 K 0, [l(| as mentioned earlier. The higher 
temperature phase has tetragonal P4/nmm symmetry 
and the lower temperature one has monoclinic P112/n, 
or equivalcntly, orthorhombic Cmma 29]. 

To optimize the structure of LaOFeSb as well as 
LaOFeP and LaOFeAs, we choose a \[2 x \/2 x 2 su- 
percell of the original tetragonal (lxlxl) primitive 
unit cell for the starting structure and fully relax the 
structure including the cell parameters. We choose the 
doubled unit cell along the c-axis to examine the exper- 
imentally observed spin configuration in LaOFeAs that 
the spin direction is opposite between two adjacent Fe 
layers along the c-axis [9(]. With this choice of the su- 
percell, we perform both the NM calculation and spin- 
polarized one with the well-known SDW spin configura- 
tion as represented in Fig. 1. 
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The optimized cell parameters and atomic coordinates 
of LaOFeSb are listed in Table I along with those for 
LaOFeP and LaOFeAs for comparison. For all of the 
compounds, our NM calculations result in tetragonal 
structures with a = b, and the lattice parameters of 
LaOFeSb are slightly larger than those of LaOFeAs as ex- 
pected. The NM-phase tetragonal structure is not stable 
with respect to AFM spin orderings. Our spin-polarized 
calculation for the SDW phase leads to an orthorhombic 
unit cell with a ^ b, and again LaOFeSb is found to have 
larger lattice parameters than LaOFeAs. Our result for 
LaOFeAs is consistent with the previous calculation [29| 
which shows that, in the SDW spin configuration, paral- 
lel spins tend to get closer and opposite spins move apart 
by structural distortion with slightly increased 7 angle 
(between a- and b-axis in the original lxl unit cell) 
from 90 degrees to lift the magnetic frustration. The 7 
angles from our structure optimization are found to be 
90.7, 90.9, and 91.1 degrees, for Pn=P, As, and Sb, re- 
spectively. The angle of 90.9 degrees for LaOFeAs is in 
reasonable agreement with the measured value of 90.3 de- 
grees supporting the validity of our predicted value 
for LaOFeSb. Moreover, the increasing 7 angle from P 
to Sb is also consistent with the fact that the calculated 
Fe spin moment increases monotonically from P to Sb, 
as presented below. 

The AFM spin configuration (Fig. 1) in the optimized 
orthorhombic unit cell is more stable than the NM phase 
in the optimized tetragonal unit cell by 153, 354, and 
706 meV per formula unit for phosphide, arsenide, and 
antimonidc, respectively. The local magnetic moment 
on a Fe atom is also increasing as 2.42, 2.83, and 3.13 
He- Therefore, the magnetism is greatly enhanced in the 
LaOFeSb compared with LaOFeAs and LaOFeP. In the 
meanwhile, no ferromagnetic (FM) solution is found in 
our calculation. 

Figure 2 shows the band structures calculated in the 
NM phase for the series of oxypnictides. Overall features 
are similar among Pn species. Complicated Fe d-derived 
bands exist near Ef, electron and hole pockets around M 
and r points, respectively, and dispersionless bands be- 
tween r and Z points, indicating the two-dimensionality 
of the layered structure. However, band widths and de- 
tailed energy positions of Fe d-derived bands are different 
among pnictides. 

As shown in Fig. 2, the band dispersion is considerably 
different among the compounds. The P compound has 
largest dispersion and band width, and then follow the 
As and Sb compounds, as is easily seen, for examples, for 
the electron bands at M point near Ef. The decreasing 
dispersion from P to Sb compounds are easily understood 
by the increasing Pn atomic size and distance between Fe 
and neighboring Pn atoms, which results in reduced or- 
bital overlap between them and therefore the localization 
of Fe d states. The localization of Fe d states close to the 
atomic orbital limit enhances the intra-atomic exchange 
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FIG. 2: (color online). Calculated band structures of a) 
LaOFeP, b) LaOFeAs, and c) LaOFeSb for NM phase. In (d), 
DOS are shown in red, green, and black colors for LaOFeP, 
LaOFeAs, and LaOFeSb, respectively. 



coupling, or, the Hund's rule coupling, so that the Fe lo- 
cal spin moment and the stability of magnetic solution 
are expected to increase from P to Sb. Our results pre- 
sented above are consistent with this expectation. 

Detailed energy positions of Fe d-derived bands are dif- 
ferent among pnictides (Fig. 2), and this affects the num- 
ber of hole pockets around the Y point. Energy splittings 
between two d z 2 -derived bands (which have finite disper- 
sion along the T-Z line), and between d xy and d yz + d zx 
derived bands (which are dispersionless along the T-Z 
line) decrease from P to Sb. As a result, there are three 
bands crossing the Ef near the T point in the case of 
LaOFeSb while there are only two bands in the cases of 
LaOFeP and LaOFeAs. The 3-D hole pocket centered 
around the Z point, which is present in LaOFeP and 
LaOFeAs, is absent in LaOFeSb. This will be discussed 
in detail later with Fig. 3. 

Figure 2(d) shows the density of states (DOS) of 
LaOFeP, LaOFeAs, and LaOFeSb in the NM phase. The 
value of DOS at Ef is 2.9 states/eV per formula unit 
in LaOFeSb, which is much greater than 1.7 states/eV 
in LaOFeAs. This implies that SDW formation and su- 
perconductivity may occur more effectively in LaOFeSb 
than LaOFeAs, since large DOS at Ef often results in 
stronger instability of the Fermi surface. In the case of 
LaOFeP, DOS at Ef is as large as 4.5 states/eV due 
to a very narrow peak right at Ep originating from van 
Hove singularities which are away from the T or M point, 
so substantial part of DOS may not contribute strongly 
to the superconductivity within the spin-fluctuation sce- 
nario. In the case of the SDW phase, the electronic struc- 
ture near Ef is drastically changed from the NM phase. 
The DOS at Ef is greatly reduced to 0.23 states/eV for 
the SDW phase due to the SDW gap formation. The 
SDW phase is still metallic, however, with a finite DOS 
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FIG. 3: (color online). Calculated Fermi surfaces in the NM 
phase: (a) LaOFeSb using the primitive unit cell, and (b) 
LaOFeSb, (c) LaOFeAs, and (d) LaOFeP using the \/2 x \/2 
supercell. Hole pockets are represented in blue, and electron 
pockets are in red. In (b)-(d), the dashed lines represent the 
Brillouin zone of the V2 x \/2 supercell. 



an enhanced nesting between the surfaces, together with 
the enhanced DOS at Ef, suggest that this material may 
have stronger pairing mechanism within the spin fluctu- 
ation scenario which can, in turn, result in superconduc- 
tivity with higher T c . 

In summary, we have investigated the basic physical 
properties of a hypothetical iron-based oxypnictide com- 
pound LaOFeSb by first-principles calculations. It is pre- 
dicted that the compound favors to be in the tetrag- 
onal crystal structure with a — b in the NM state 
while orthorhombic structure with a ^ b is the lowest 
energy structure in the stripe-type AFM-SDW phase. 
The SDW phase becomes more stable, having the local 
spin moments on Fe atoms and the total energy differ- 
ence between the SDW and NM phase larger than those 
in LaOFeAs. In the NM phase, the band dispersion 
around Ef is substantially reduced, increasing DOS at 
Ef greatly. The calculated Fermi surface in the NM- 
phase LaOFeSb exhibits enhanced nesting feature be- 
tween the hole and electron surfaces with the nesting 
vector (7r,7r) for the SDW formation. 

This work was supported by the KRF (KRF-2007-314- 
C00075) and by the KOSEF Grant No. R01-2007-000- 
20922-0. Computational resources have been provided by 
KISTI Supercomputing Center (KSC-2008-S02-0004). 



value at Ef, as the nesting is not perfect as in the case 
of LaOFeAs system . 

The Fermi surface of LaOFeSb is shown in Fig. 3 
on the k z — plane. The Fermi surfaces of LaOFeAs 
and LaOFeP are also shown for comparison. Similarly 
to the LaOFeAs and LaOFeP cases, there are hole pock- 
ets around the F point and electron pockets near the M 
point in LaOFeSb. In this case, however, there is ad- 
ditional small hole Fermi surface centered at F which is 
absent in LaOFeAs and LaOFeP, as a result of different 
Fe-d band distributions near Ef (see Fig. 2). There- 
fore, LaOFeSb has three hole pockets near the F point 
and two electron pockets around the M point. All Fermi 
surfaces are cylindrical. The small hole surface at F has 
different orbital character from the 3-D hole pocket at 
Z in LaOFeAs and LaOFeP. The nesting between elec- 
tron and hole Fermi surfaces defined by the (7r,7r) vector 
in the reciprocal space is also manifest in this case. In 
order to make comparison straightforward, Figs. 3(b)- 
(d) show the original Fermi surfaces and those shifted by 
the (±7r,±7r) vectors. As shown in Fig. 3(b), the shapes 
of Fermi surfaces in LaOFeSb are quite close to circles 
compared with those in LaOFeAs and LaOFeP, and the 
radii of circular surfaces match closely between the hole 
and electron surfaces. Furthermore, the electron surfaces 
coincide with each other almost perfectly, and so do the 
two hole surfaces at F. These features which indicate 
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